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Abstract 

We study the adiabatic density perturbation in the oscillating inflation, proposed by Damour and Mukhanov. The 
recent study of the cosmological perturbation during reheating shows that the adiabatic fluctuation behaves like as 
the perfect fluid and no significant amplification occurs on super-horizon scales. In the oscillating inflation, however, 
the accelerated expansion takes place during the oscillating stage and there might be a possibility that the para- 
ge ■ metric amplification on small scales affects the adiabatic long-wavelength perturbation. We analytically show that 
^\ \ the density perturbation neglecting the metric perturbation can be amplified by the parametric resonance and the 
i-H ■ instability band becomes very broad during the oscillating inflation. We examined this issue by solving the evolution 
^ ' equation for perturbation numerically. We found that the parametric resonance is strongly suppressed for the long 
wave modes comparable to the Hubble horizon. The result indicates that the metric perturbation plays a crucial role 
for the evolution of scalar field perturbation. Therefore, in the single field case, there would be no significant imprint 
of the oscillating inflation on the primordial spectrum of the adiabatic perturbation. However, it could be expected 
that the oscillating inflation in the multi-field system gives the enormous amplification on large scales, which may 
lead to the production of the primordial black holes. 
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' I. INTRODUCTION 

<n : 

Dynamics of the coherently oscillating scalar field play an important role in the early stage of the universe and has 
been studied by many authors. Specifically, the significance of parametric resonance is found in the recent investigation 
of the reheating process (!]]. Due to the coherent oscillation of the background inflaton field, the light boson field or 
the fluctuation of the inflaton itself is amplified through the non-linear interaction or the self- interaction ||||. The 
efficiency of the parametric resonance could have the important implication for GUT scale baryogenesis . 

From the viewpoint of the large scale structure in the early universe, several authors investigated the influence of 
parametric resonance on the primordial density perturbation during the coherently oscillating stage after inflation. 
In the framework of cosmological perturbation, the single field model is studied in Ref. |5|||. The analysis is also 
extended to the two-field model by Taruya and Nambu [Q] . These analyses are mainly focused on the long wavelength 
perturbation, which significantly differs from the one neglecting the metric perturbation. Nambu and Taruya j^] 
investigate the Mukhanov's gauge-invariant variable and found that the evolution equation for perturbation can be 
reduced to the Mathieu type equation. Although the Mathieu equation itself has the exponential instability, the 
instability of the perturbation is very different from this and the final result is fully consistent with Kodama and 
Hamazaki 0). The important conclusion of these analyses is that the density perturbation on super- horizon scales 
behaves like as the perfect fluid unless the iso-curvature perturbation becomes dominant. Therefore, as far as the 
adiabatic perturbation is concerned, the primordial power spectrum on large scales does not suffer from the significant 
amplification by the parametric resonance. However, there remains a possibility that the parametric resonance can 
appear well inside the Hubble horizon. For example, consider the massless self-interacting inflaton. Neglecting the 
metric perturbation, it is known that the fluctuation of the inflaton field has the instability mode whose wavelength 
is much smaller than the Hubble horizon . 

The aim of the present paper is to clarify the perturbaton on sub-horizon scales taking into account the gravitational 
perturbation and understand the cosmological implication of the parametric amplification to the large scale structure. 
As usual, the universe expands decelerately during the coherently oscillating stage and the parametric resonance inside 
the Hubble horizon cannot affect the density perturbation on super-horizon scales. However, Damour and Mukhanov 
recently considered the model with the non-convex type potential S. In this model, taking the time average, the 
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potential energy of the inflaton becomes large compared to the kinetic energy and the accelerated expansion can take 
place in the oscillating stage. Liddle and Mazumdar called it oscillating inflation and confirmed this fact numerically 
p0| | . As was suggested by Damour and Mukhanov, it could be expected that the very broad parametric resonance 
occurs and the quantum fluctuation is enormously amplified. This implies that the adiabatic metric perturbation with 
the large amplitude might be produced on super-horizon scales in the course of the accelerated expansion, which can 
lead to the different conclusion from the previous results H[|. Hence, there exists a possibility that the cosmological 
objects such as the primordial black holes may be formed by the large amplitude of the metric perturbation. 

In this paper, to explore this possibility, we study the cosmological perturbation in the oscillating inflation and 
investigate the efficiency of the parametric resonance. We describe the model in Sec|n]. The parametric resonance 



of density perturbation is discussed in Sec. HI. We show that the efficiency of the parametric amplification on the 
spectrum of the curvature perturbation depends on the energy scale of the oscillating inflation and the growth factor. 
Section IV is devoted to the analysis of the growth factor. Neglecting the metric perturbation, we analytically find 
that the instability band of the parametric resonance is very broad and the growth factor may become significantly 
large. We then proceed the numerical calculation of cosmological perturbation without any approximation. The 
numerical analysis show that the perturbation can experience the parametric amplification, however, the amplification 
is strongly suppressed when the wavelength of the fluctuation approaches the Hubble horizon scales. The results are 
briefly summarized and the conclusion in the single field case is described in the last section [v|. We also discuss the 
oscillating inflation in the multi-field system. It could be expected that the significant amplification of the density 
perturbation can appear on large scales and may lead to the production of the cosmological black holes. 



II. OSCILLATING INFLATION MODEL 

Let us consider the minimally coupled scalar field in the flat FRW universe. The homogeneous background 
equations become 

pi 

4> + 3H$ + V. 4) = 0, (2) 

where H is the Hubble parameter and M p i is the Planck mass. The evolution of the scalar field depends on the shape 
of the potential V ((/>). We shall consider the model of oscillating inflation proposed by Damour and Mukhanov ||, 
in which the potential is given by 



V{4>) = - 



q/2 



(3) 



Note that the potential (||) becomes logarithmic, V = {A/2) log ((f) 2 /<p 2 + 1) in the limit q —> 0. 

When the effective mass of the inflaton is larger than the Hubble parameter, i.e, {V^l ^> H 2 , the inflaton field 
shows oscillatory behavior. For the large amplitude 4> 3> </> c , taking the time average during a period of oscillation 
yields the following virial theorem: 

(^)^q(Vm. (4) 

Using this relation, the inflaton field can be regarded as the perfect fluid matter with the equation of state P = 
(lose — l)p ; where j osc is the effective adiabatic index given by 

l0SC = — — y (5 

Thus the time dependences of the scale factor, the amplitude of the inflaton field and inflaton potential are obtained 
§ & 

a(t) cx ^+ 2 )/ 3 ?, 

4>{t) oc a - e/< - q+2) cx r 2 /\ (6) 
V{4>) cx a - 6 9/("+2) K t -» i 

where d> is the amplitude of the scalar field evaluated at the time that = holds. 

Eq.(ph implies that if the potential index q is smaller than unity, the universe can experience the accelerated 
expansion even in the oscillating phase. When the scalar field approaches to the critical value ~ </> c , the oscillating 
inflation ends. Hereafter, we shall restrict our attention on the inflationary phase with the potential index q < 1. 



2 



III. COSMOLOGICAL PERTURBATION 



We now investigate the cosmological perturbation in the oscillating inflation. We shall ignore the interaction with 
the other scalar fields which leads to the inflaton decay in the reheating process In Ref. the useful and 

well-behaved perturbed quantity during the oscillating stage is found : 



Q = 



(7) 



where S(j) and 1Z are the perturbation of the scalar field and the spatial curvature, respectively The quantity Q 
is the gauge-invariant variable introduced by Mukhanov and the evolution equation is simply given by 



Q + 3HQ 



V, 



Q = 0. 



(8) 



If we neglect the third term in the bracket, the above equation is just reduced to the equation for perturbation 84> 
ignoring the metric perturbation. We thus understand that the effect of gravitational interaction to the scalar field 
perturbation is encoded in the third term in the bracket. 

In equation (@), there exists the exact solution in the long- wavelength limit Hj. For k — > 0, we have 



Q = d 



H 



H 



dt H 2 



(9) 



where c\ and C2 are the integration constants. The solution proportional to the coefficient c\ is named as the growing 
mode and the one proportional to C2 as the decaying mode. [] For the growing mode, the curvature perturbation on 
comoving slice 1Z C defined by 



H 



(10) 



remains constant in time. The curvature perturbation is also referred to as the Bardeen parameter, which is simply 
related to the CMB temperature fluctuation AT/T observed by the Cosmic Background Explorer satellite [^9|. In this 
paper, matching the amplitude of quantum fluctuation Q produced during the oscillating inflation inside the horizon 
with that of the long wavelength growing mode, we will evaluate the power spectrum of the curvature perturbation 
V n {k). 

We focus on the evolution of the quantity Q on sub- horizon scales. In equation (p|), the time dependence of the 
terms in the bracket is estimated by using the results (0): 



V tH> (0) 



H 



-6(2=1} 



(11) 
(12) 



which are valid for ^ > </i c . During the oscillating inflation, the second term in the bracket becomes dominant 
compared to the term ([l2]) . The dominant contribution of the oscillating term ( fj"l| ) differs from the situation considered 
in Ref. || Q. We then introduce the following variables: 



drj = g(t)dt ; g(t) 



q/2 



r 3 (f+t). 



Ignoring the term (12), the evolution equation m) is rewritten as the normal form 



t It seems that the singular behavior appears in the decaying mode around the zero points of the time derivative of the scalar 
field. However, the zero points of the denominator in the integral can be canceled out by the zero points of <j> in the numerator. 
Therefore the solution is regular and the amplitude of the solution decreases in time. For more rigorous proof of the regularity 
and the explicit calculation of the long- wavelength perturbation, see [jq] and M . 
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-) +v«W 



= 0, 



(13) 



where is defined by <f>(t) = <f> 4>[rf) and satisfies the equation of motion 

km + V$) = °' ( 14 ) 

with the initial condition <j>{r)i) = 1. The quantities V ^ and V ^ are obtained by differentiating the effective potential 
V with respect to cj): 



2\ 9/2 



(15) 



Since we have V ~ ^VoS <^ oscillates rapidly during the oscillating inflation. Treating the time dependence of the 
term V adiabatically, the period of oscillation 4> is estimated as 



T 




(16) 



'2(g-i - V) 



For the range 0.01 < q < 1, we have T ~ 5. 

According to the Floquet theorem, presence of the oscillating term i n (|13| ) implies that there exists the following 
solution for the equation (Oh. Since V it is given by 77 = T/2, we have plfl 



Q(r) = e» T P(r) ; P(r + tt) = P(r), 



(17) 



where we define the new time parameter t = (2tt /T)r). The periodic function P(t) is bounded and has a finite 
amplitude. The characteristic exponent fj, takes either a imaginary or a real number, which depends on the shape of 
V 21 and the parameter k/ga. The parametric resonance means that the characteristic exponent has a real number. 

We can see that the parametric resonance can affect the curvature perturbation 1Z C produced during the oscillating 
inflation as follows. As usual, in the slow-rolling inflation, the amplitude of quantum fluctuation is given by Q ~ 
\5<ft\ ~ H/2tt, which can be deduced from the instability of the massless scalar field in de Sitter space jl5) . On the other 
hand, the amplification of the quantum fluctuations can appear due to the parametric resonance in the oscillating 
inflation. We have 



(18) 



where subscripts (0) and (*) denotes the quantities evaluated at the beginning of oscillating inflation and at the 
time when the wavelength of the quantum fluctuation exceeds the Hubble horizon scale (k = a* if*), respectively. 
Here, At denotes the time interval during which the quantum fluctuation is amplified by the parametric resonance. 
For the modes inside the horizon k/aH — * 00, the mode function of the equation (H) initially becomes etoQo ~~ * 
exp (i J dt/a)/V2k, corresponding to the vacuum state in the Minkowski spacetime | fl6| . Matching the amplitude of 
the short wavelength fluctuation with that of the long wavelength solution (||), the curvature perturbation (J^) is 
roughly estimated : 




(4-g)/(2+<7) 



Qo 



(19) 



which is evaluated at the horizon-crossing time. Therefore, the power spectrum of the curvature perturbation Vn(k) 
evaluated by taking the ensemble average (• • -) en s becomes 



2tt 



(A-q)/(2+q) 



V, 



1/4- 



M, 



(20) 



pi 



In the derivation of the expression (|2C| ), we have naively replaced the ensemble average (4> 2 ) en s with the time average 
{(j) 2 ) and used the virial theorem (|4|). Eq.(|2(]) shows that the appearance of exponential factor e MAr may give the 
important contribution to the spectrum of curvature perturbation, although the amplitude of V-ji is always suppressed 
by the energy scale of the oscillating inflation V^ 4 = V /1 / 4 (0») |l(| |l~9|]. 
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IV. PARAMETRIC AMPLIFICATION DURING OSCILLATING INFLATION 



In the previous section, we have seen that the effect of parametric resonance can lead to the amplification of 
the curvature perturbation, which may gives V-ji ~ 1. In this section, to explore the efficiency of the parametric 
amplification, we analytically and numerically investigate the growth factor e AlAr . 



A. Analytic estimation 



Let us evaluate the characteristic exponent fj, from the reduced equation (13). Using the solution ([l7|) obtained 
from the Floquct theorem, we have the following formula H]: 



fj, = - log 



^ + Vf 2 - 1 



F=l + 2Qi(r = |)Q 2 (r=|; 



(21) 



where the prime denotes the differentiation with respect to the time r. Qi and Q2 are the solutions satisfying the 
initial conditions Q\ = 1, Q\ = and Q2 = 0, Q' 2 = 1 at r = 0, respectively. The parametric resonance appears in 
the case F > 1. 

Although it is difficult to obtain the solutions Q\ and Q2 for our complicated potential V11, we can get the 
approximate expression of the characteristic exponent. The effective mass of the oscillating innaton field becomes 
negative at the shallow wings of the potential and positive at the minimum. Thus we replace the periodic function 
V 77 with the step-wise function given by 




(22) 



where we choose 



2tt 

n = Y 



d<p 



2{tr 1 -v) 



1-9 



At <fi — 0i, we have — 0. The approximation ( p2| ) enables us to obtain the analytic solutions Qi and Q 2 - Then 

the quantity F is evaluated as follows: 



2 \p a J 2 2 



(23) 



where 



T 
2^ 



1/2 



T 



1 1/2 



-) +v«(* = i) 



Since V ^{4> = 1) < during the oscillating inflation, the variable /? becomes imaginary for the long- wave modes 

(k/a) 2 < |K</>0(^)l an( i we obtain F > 1. This can be deduced from the negative coupling instability ||. Important 
point is that we could also have F > 1 for the short-wave mode (k/a) 2 > |U00(0)|. In the latter case, the band 
structure of the characteristic exponent /1 appears. 

In Fig.l, we plot the characteristic exponent /1 as a function of k/ga evaluated from ( ^l| ) and (p3|). The parameters 
are chosen as C = 10~ 6 M p / and q = 0.1. As for the quantity 4>, we treat it as constant and set </> = q M p i/\/6, 
corresponding to the initial value of the oscillating inflation |l0| . The circles show the characteristic exponent obtained 
from the substitution of the numerical solutions of equations ([l3]) and ( |l4| ) into the formula (|2l"|). We see that the 
approximation predicts the numerical result reasonably well. The crucial observation is that the instability band is 
very broad and the characteristic exponent can be of the order of unity. 

To evaluate the growth factor e MA , we further need the period Ar. It is given by 



Ar 



Y I ''''Hi) 



T2 q{q 



r(| + |) ( m, 
r(|) 



P i 



(24) 
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where the integral is evaluated from the time at which the oscillating inflation starts to the horizon-crossing time. 
Eq.(^4|) is valid for the case q > 0. Rewriting the above equation with the relation k = a*H*, the wave number 
dependence of the growth factor e MAr shows that the characteristic peak appears on the spectrum of curvature 
perturbation (|2p|). Because of the rapid oscillation of the inflaton field, the period At effectively becomes very large. 
For q = 0.1, the typical values 0* = 10" 4 , 1CT 5 and 1CT 6 M pt give At = 1367, 13696 and 136991, respectively. 
Therefore, even if we get the rather small value of the characteristic exponent [i = 0.01, for example, the contribution 
of the growth factor to the curvature perturbation becomes significantly large. For 0* = 10 _5 M p j and q = 0.1, we 
have e^ AT = 10 182 ! 



B. Numerical result 



The analytic estimation in the previous subsection rather overestimates the growth factor. For more rigorous 
evaluation, we must take into account the time dependence of <p an d k/ga. Furthermore, we should remember that 
the previous result comes from the analysis of the reduced equation (|l3|) . This also leads to the overestimation of the 
growth factor. However, there still exists a possibility of the parametric amplification of the curvature perturbation 
spectrum. To clarify the influence of parametric resonance, we numerically solve the equations (0), ^ and (||). 

Fig. 2 shows the growth factor e^ AT " for the index q — 0.1, i.e, the ratio of the amplitude Q e evaluated at the end 
of oscillating inflation to the amplitude Qo at the beginning of inflation by varying the parameter <j) c . For each value 
(j) c , we start to calculate the background evolution by setting the value cf>o — qM p i/^/& corresponding to the final 
value of the slow-rolling inflation. As for the fluctuations, the initial conditions Q — are chosen. Although we have 
checked the evolution for the various initial conditions Qo cx cos (k J dt/a), corresponding to the vacuum state in the 
short wavelength limit, the maximal amplitude can be obtained from the initial condition Q' = 0. Fig. 2 shows that 
the amplification for the fluctuations occurs by the effect of parametric resonance. It is obvious that the maximal 
amplitude depends on the parameter 4> c . For the smaller value of 4>c, the amplitude of the fluctuation becomes larger, 
which comes from the fact that the period of inflation represented by the e-folding number N e = log e (a e /ao) depends 
on the critical value cf> c |]|l(|. When the oscillating inflation takes place for a longer time, the effect of parametric 
resonance works out more efficiently and the fluctuations are significantly amplified. In Fig. 2, the three different lines 
are depicted. The solid line is the fluctuation with the wavelength ten times smaller than the radius of the Hubble 
horizon at the end of oscillating inflation, i.e, k = 10a e H e . The dashed line has the wavelength k = 2a e H e and 
the dotted line corresponds to the fluctuation whose wavelength just reaches at the Hubble radius after oscillating 
inflation. It is remarkable that the modes inside the Hubble horizon are enormously amplified, while the amplification 
of the long- wave mode comparable to the Hubble horizon is strongly suppressed. Though the dotted line apparently 
takes the growth factor ~ 10 7 , it can be recognized that the parametric resonance becomes ineffective for the long-wave 
modes. 

To see the efficiency of the parametric amplification explicitly, we plot in Fig. 3 the time evolution of the fluctuations 
using the quantity Q from the beginning to the end of oscillating inflation. We set the parameters q = 0.1 and 
4> c = lO~ 5 AI p i. The horizontal axis represents the cosmic time normalized by the factor \[AjM v i. Each figure 
shows the time evolution with the different wavelength under the same initial conditions Qo = 1.0 and Q' — : (a) 
k = 10a e H e ;(b) k = 2a e H e ; (c) k = a e H e . These figures clearly reveal that the growth of the perturbation is much 
sensitive to the wavelength of the fluctuations. The parametric amplification works efficiently inside the horizon, 
however, the amplitude Q becomes constant for the wavelength near the Hub ble horizon. 

We can understand these behavior as follows. As was described in Sec. Ill, there exists the exact solution for the 
variable Q in the limit k — > 0. Ignoring the decaying mode proportional to the coefficient C2 in (^), the amplitude 
of the long wavelength solution becomes nearly constant, which can be deduced from (||). This behavior can also be 
obtained in the case of k ^ 0, when the term (k/a) 2 in equation (||) is negligible compared to the terms and 
M ~ 2 (V/H)' , which gives the condition k/a <C H ||. In the previous subsection, we analyzed the growth factor e AtAr 
for the short-wave mode by dropping the term M~ 2 (V/ H)' . Fig. 3 indicates that the term induced by the gravitational 
perturbation plays a crucial role even for the modes comparable to the Hubble horizon. This is fully consistent with 
the result of the paper |^] and Therefore, the spectrum of the curvature perturbation produced during oscillating 
inflation on super-horizon scales has the rather small amplitude compared with the analytic prediction. 



V. CONCLUDING REMARKS 



We have analyzed the cosmological perturbation in the oscillating inflation and explored a possibility of parametric 
amplification for the curvature perturbation. Neglecting the metric perturbation, the analytic estimation shows 
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that the curvature perturbation could be amplified by the broad band parametric resonance. We then numerically 
examined this issue by solving the evolution equation for perturbation without any approximation. The enormous 
amplification of the perturbation can appear, however, we found that the presence of the the metric perturbation 
strongly suppresses the parametric amplification on large scales comparable to the Hubble horizon. 

Now we discuss the cosmological implication of the parametric resonance to the density perturbation. According 
to Ref. Jl8],|l9| , the physical length of the fluctuation produced during the oscillating inflation is 

X P hy < 10 1+OA3a * cm, (25) 

where X p h y denotes the present wavelength. The quantity a* is given by 

a.= ^+log e p° 16 1 G 4 eV V (26) 



v}' 4 /' 

iV* denotes the e-folding number of the oscillating inflation evaluated at the Hubble crossing time. On the other hand, 
if the amplitude 1Z C produced during the inflation becomes of the order of unity, the fluctuations can experience the 
gravitational collapse. Assuming that the primordial black holes are formed at the horizon re-entry time during the 
radiation epoch, we can evaluate the typical mass of a black hole using the equation (^|j). We obtain 

Mbh < f0 1+a86Q * g, (27) 

which is evaluated at the formation epoch. 

As is shown by Liddle and Mazumdar |T(|, a period of oscillating inflation is short. This can be checked in our 
numerical calculation. For the critical value <fi c = f Q~ 6 M p i, we obtain the small value of the e-folding number N e = 3.5. 
Hence, the expressions ( p5| ) and (|27j) show that the energy scale of the oscillating inflation should be rather low to 
have the cosmologically interesting scale for the back hole mass and the physical size of the fluctuations. Asaka, 
Kawasaki and Yamaguchi |2(J investigated the inflation model to resolve the cosmological moduli problem, in which 

the oscillating inflation occurs at the energy scale Ve^ 4 — 1 OGeV, instead of the thermal inflation. In this energy 
scale, if we have the e-folding number N e ~ 3, it is possible to get Mbh ~ LM©j as a candidate for explaining the 
recent observation of the massive compact halo objects, although the physical length is still small, X p h y ~ 1 pc. 

However, the amplitude of the vacuum fluctuation becomes extremely small for the inflation at low energy scale. 
We saw in Sec. HI that the spectrum of curvature perturbation Vn given by ( f20| ) contains the suppression factor 



(V* 4 /Afpj) 2 , which takes ~ I0~ 36 in the above case. Furthermore, our results show that the gravitational interaction 
plays a crucial role for the parametric amplification of the curvature perturbation. To obtain the relevant growth 
factor e^^ T to compensate the suppression factor, a long period of oscillating inflation is necessary, which requires 
the very fine-tuning of the model parameter </> c . Therefore, in the single field case, we conclude that the influence of 
the parametric resonance on the cosmological perturbation would not be imprinted on the universe observed as long 
as the adiabatic mode of the perturbation is dominant. 

The conclusion might be changed if we consider the multi-field system. Consider the hybrid-type inflation which 
is a more realistic model motivated by the particle physics Jf9| . In this scenario, we can imagine that the oscillating 
inflation occurs subsequent to the slow-rolling inflation and it is followed by a secondary inflation driven by another 
scalar field. If we choose the appropriate duration of the secondary inflation, the oscillating inflationary phase could 
have the relevant length scales \ p h y which can affect the large scale structure. In addition, the hybrid inflation would 

make a modification of the spectrum Vn- In equation j20|), we should replace the term (V^ 4 /M p i) 2 with V*^/M p i 
|P7[ . This relaxes the fine-tuning of the model parameter </> c to give the cosmologically interesting scale of the black 
hole mass Mbh- The effect of oscillating inflation might be imprinted on the primordial density perturbation in 
such inflationary scenario. To investigate the spectrum Vn of the hybrid inflation, we must study the cosmological 
erturbations with multi-field system, in which the contribution of iso-curvature perturbation cannot be neglected 
The non-linear interaction of the oscillating inflaton with the secondary inflation-driven field would lead to the 
significant amplification of the iso-curvature perturbation by the parametric resonance. Although the influence of 
iso-curvature mode on the gauge-invariant perturbation 1Z C has not been understood completely, the result in this 
paper would shed light on the evaluation of the primordial power spectra. We will report the analysis of the multi-field 
system in a separate publication. 
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FIGURE CAPTION 



Figl Characterisctic exponent \i for the parameters q — 0.1 and <fi c — 10~ 6 M p i. We plot the exponent /i as a function 
of the Fourier mode k/ga. The solid line shows the approximation using the result (|23|). The circles are the 
characteristic exponent evaluated by substituting the numerical solutions of eqs.(O) and (|14|) into the formula 



one Qo at the beginning of inflation. We plot the ratio of the amplitude with the same initial condition Q' Q = 
by varying the model parameter </> c . We set the potential index q = 0.1. The solid line corresponds to the 
fluctuation with the mode k = l0a e H e , i.e, the wavelength is ten times smaller than the Hubble horizon size at 
the end of oscillating inflation. The dashed line represents the amplitude for the modes k = 2a e H e . The dotted 
line is the ratio for the fluctuation whose wavelength reaches at the Hubble horizon size just after the oscillating 
inflation ends. 

Fig3 Evolution of Q in the case of the parameters q = 0.1 and 4> c — 10 _5 M P ;. For each figure, we start to calculate 
the background eqs . (Jl|) and (||) from the value <f> = q/VQ after slow-rolling regime. The initial conditions for 
the fluctuation Q are set by Q — 1.0, Q' — and eq.(pt) is solved numerically. The horizontal axis denotes the 
cosmic time normalized by y/A/M p i: (a) The fluctuation Q for the mode k = 10a e H e ; (b) The fluctuation Q 
for the mode k = 2a e H e ; (c) The gauge-invariant quantity Q whose wavelength reaches at the Hubble horizon 
size just after the oscillating inflation ends, i.e, k = a e H e . 
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Fig.2 
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Fig.3c 
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